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Chapter 1 


INTRODUCTION 


The release of volatile organic compounds (VOC’s) into 
the atmosphere has come under fire by the Environmental 
Protection Agency (BPA) and state regulatory agencies in the 
past several years. These releases come from a wide variety 
of sources from gasoline filling stations to industrial and 
military operations. There are also countless natural 
sources of VOC’s including volcanos and forest fires. 
Environmental groups have had a great influence on 
government agencies resulting in the enactment of numerous 
statutes and regulations to reduce these releases over an 
approximate five (5) year period. 

Enforcement of environmental regulations is carried out 
by the individual states which impose their own specific 
regulations which in some cases are even more stringent than 
those of the federal government. A case in point is the 
State of California where regulations are resulting in rapid 
actions on the part of corporations to comply or face harsh 
penalties. These actions include both a switch to 
alternative manufacturing materials and employment of new 
pollution control devices. 

A new pollution control technology currently under 
development utilizes ultraviolet (UV) light to break down 


the volatile organic compounds (VOC’s) contained in an 





2 
airstream. One pollution control equipment manufacturer has 
already employed this method in an installation in 
California (see Figure (1)). The process involves injecting 
ozone into a mixture of gaseous VOC’s which have been 
radiated by UV light. In the presence of oxidants and UV 
light, the air/VOC mixture undergoes a series of photo- 
dissociation and radical oxidation reactions which destroy 
the VOC’s. Despite this pollution control equipment’s 
success in controlling a wide range of airborne VOC’s, the 
extent of VOC reduction due to photolysis in the photolytic 
reactor section of the system is unknown. Our lack of 
information on the effectiveness of this photolytic reactor 
as shown in Figure (2) is due to proprietary information 
retained by the manufacturer, insufficient test data, and 
lack of an adequate model to predict its performance. 

There have been some analytical and experimental 
Studies reported in the literature on photochemical reactors 
(l1e., Y. Harano and T. Matsura (1972)). However, these 
studies are generally directed towards simplified geometries 
which can be modeled analytically (ie., axial, elliptical, 
and columnar reactors) and which assume well-mixed 
Somait tons, low elagne abseorpcion, uniform flow, and 
negligible diffusion. Experimentation has been directed 
towards determining the deviation between these highly 
Simplified models and actual reactor performance in order to 


determine the effects of factors omitted. Few if any 
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TAES ULTRAVIOLET-OXIDATION CONTROL SYSTEMS have been designated as Best 
Available Control Technology (BACT) for Volatile Organic Compounds (VOC) and Reactive Organic 
Compounds (ROC). 


The UV-Oxidation processes utilize a proprietary, on-site, self-regenerative system designed to meet 
the current EPA and most stringent local Air Quality District rules for capture and destruct 
efficiencies. 
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Figure (1) - VOC Pollution Control Device Schematic 


(Terra-Aqua Environmental Systems (1992) ) 
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Figure (2) - VOC Pollution Control Device Photolytic Reactor 





studies can be found in the literature which attempt to 
quantitatively model photochemical reactors of complex 
geometries and flow fields using numerical methods. This 
report provides a numerical model which can be used to 
quantitatively describe the photo-dissociation/radical 
Graligatalon GCE VOC san the photolytic reactor of Figure (2). 

A secondary but equally important objective for this 
analysis is to develop an analytical tool which can be used 
to optimize photo-reactor geometry. This numerical model, 
once verified with experimental data (in progress at Penn 
State) and expanded to incorporate additional reactants, can 
be utilized to predict optimum geometry for the photolytic 
reactor section of a pollution control device. 

Finally, this model will assist the designer in 
determining when the use of a photolytic reactor ina 
pollution control device is appropriate. Since there are a 
multitude of VOC’s which vary widely in their ability to 
react when exposed to UV light, cost effective 
implementation of a photolytic reactor may be highly 
dependant on the specific VOC’s expected in an airstream. 
For example, the photolytic reactor may not provide any 
additional benefit in a pollution control device which 
treats specific emissions from a military or industrial 
painting operation. This numerical model may indicate when 
a significant cost savings is possible by omitting the 


photolytic reactor section and could also predict when the 





PaemetyEteereaceor Can fumction as the sole VOC reduction 


device in a pollution control system. 


1.1 Literature Review 

A review of literature on the subject of design and 
analysis of photochemical reactors shows that very little 
work has been Piasea in the area of numerical 
modeling. 

Teaeoo oe. ana DPranoftt, J. S. (1969) conducted 
experiments on light distribution within an elliptical 
photoreactor and compared these results with some closed 
form solutions. They found that the closed form solutions 
are unable to accurately predict light intensity within an 
elliptical reactor. Additional studies by Jacob, S. M. and 
Dranoff, J. S. (1970) concluded that numerical means are 
required to accurately determine the light intensity 
distribution within a well-mixed elliptical photoreactor. 
While pointing out the necessity to use numerical methods in 
the analysis of photolytic reactors, neither of the above 
Studies attempted to deal with reactors containing non- 
uniform flow fields and spatially varying reactant 
concentrations. In both studies only the case of simple 
elliptical geometry with laminar flow and well-mixed 
conditions was considered. 

Work by Ragonese, F. P. and Williams, J. A. (1971) on 


the scale-up of a laboratory photoreactor model to meet 





actual chemical production requirements outlines the 
difficulties in utilizing both dimensional analysis and 
empirical data. Although success has been achieved in this 
scale-up process, the methods used apply to simplified 
geometries with laminar flow and well-mixed conditions. 
This study points out the need for numerical modeling to 
enable a designer to scale photoreactors to meet production 
needs or to achieve sufficient pollutant reductions ina 
pollution control device. 

The literature review provided by Harano, Y., and 
Matsura, T. (1972) describes the many difficulties involved 
in the design of photochemical reactors. They summarized 
that even for the simplified geometries which have been 
studied extensively to date including cylindrical, parallel 
plate, annular and elliptical reactors, the actual 
concentration gradients, light intensity distributions and 
flow conditions greatly influence overall reaction rates 
within the photoreactor. 

All the literature reviewed on the subject of 
photolytic reactor design points to the need for numerical 
reactor modeling which incorporates both photochemical and 
fluid mechanics principles and the latest numerical solution 


techniques. 





Chapter 2 


PRELIMINARY MODELING ASSUMPTIONS 


Due to the great difficulty in modeling the entire 
reactor under investigation including all possible chemical 
reactions, this model describes only the photo- 
dissociation/radical oxidation (photolysis) of formaldehyde 
in a rectangular flow field perpendicular to the axis of a 
single UV lamp as shown in Figure (2). Formaldehyde has 
been chosen for this model since its photochemistry has been 
studied, documented and verified. Although this model 
applies only to the photolysis of formaldehyde, it is 
general enough to incorporate the chemistry of other 
airborne components as well as combinations of other 
photochemically reactive solvents in the future. In 
addition, the principles required to expand this numerical 
analysis to additional materials and/or chemical reactions 
have been included in this report. 

The following assumptions have been made to simplify 


the analysis: 


A. Constant fluid properties (ie.,density, 
temperature, diffusion coefficient). 

B. Negligible thermal effects (ie.,solution of energy 
conservation equations not required). 


C. Two-dimensional flow field. 





De Weak assorbance Gof UV light (ie., light intensity 
is independent of formaldehyde concentration). 

Hee ieOtLe-Gisse@ciation/radical oxidation of 
formaldehyde only. No other reactions between formaldehyde 
and other compounds are included except those required to 
complete the photo-dissociation/radical oxidation of 
formaldehyde. 

F. Steady, time-averaged flow field (ie, no 
oscillation of vortex shedding in the wake region). 

G. No wall losses (ie., chemical reactions which 
sometimes occur due to collisions of molecules into the 
reactor walls or the UV bulb itself). 


H. No reflection of UV light from the reactor walls. 


The above assumptions are considered to be reasonable 
ones for use in this initial numerical model. In the 
future, factors such as wall losses and reflection can be 
introduced into the numerical simulation to increase its 


accuracy . 
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Chapter 3 


MASS CONSERVATION WITHIN THE PHOTOLYTIC REACTOR 


Given the assumptions provided in Chapter 2, the 
differential equation which describes material concentration 


in the photolytic reactor can be written: 


-vove + DV2c +Q = dc/dt fate 


where, C = concentration of formaldehyde (PPM). 


Vv flow velocity (with components VurVy) (m/s). 


D = formaldehyde diffusion coefficient (m2/s) . 


(2 = generation of formaldehyde (PPM/s) . 


The first element of Equation [1], the convection term, 
involves the fluid mechanics of the reactor. It combines 
concentration changes along a given direction with the fluid 
velocity to result in a time rate of change in 
concentration. 

The second term in the above equation is the diffusion 
term. It describes the time rate of change of concentration 
entering/leaving the control volume due to concentration 
gradients at the element surface. Note that this term is 
dependant on the second spatial derivative of C multiplied 
by a diffusion constant and is relatively small with respect 


to the velocities in the flow. To achieve the highest 





ek 
accuracy possible for all possible concentration gradients 
and flow velocities, this term is included in the model. 
The third term in this equation describes concentration 
changes due to the photolytic action of the UV bulb and 
other chemical reactions. For this model the "generation 


term" simplifies to the following: 


OQ = -KC 24] 


where, K is the overall specific absorption rate constant 
(overall rate constant for the breakdown of formaldehyde by 
Bi elagnt) . 

The final term is the resultant time rate of change for 
the element which, when included yields a transient 
solution. Since the transient solution has practical value 
both in predicting the start-up concentrations within the 
reactor as well as in examining time varying inputs, the 
time dependant term has been included. In addition to the 
above benefits, a transient solution must be found in order 
to obtain convergence for the numerical solution technique 
used. This is required since it is difficult for the 
iterative numerical technique (Gauss-Seidel) to jump 
directly from zero time to an infinite solution (steady- 
state solution) without some intermediate steps in time. 

For the case of a multi-component flow field which 


includes chemical reactions as well as photolytic reactions, 





1s? 


Equations [1] and [2] can be written: 


oven bay-c, + 2 = dc,/dt boy 


Q. im eee crmMer redactions (see Appendix B) [4] 


Analysis of the multi-component case requires the solution 
of n simultaneous Equations ([{3] and [4]) corresponding to 
the total number of reactants. Note that Appendix B 
contains numerous additional equations involving chemical 
and photochemical reactions between many possible airborne 
constituents. In order to simplify this analysis only the 
photolysis of formaldehyde will be considered (ie., no 
reactions between formaldehyde and other airborne 
constituents will be considered). 

In the next chapter the determination of velocities 
required in the convection term of Equation [1] is 


discussed. 





is 
Chapter 4 


FLUID MECHANICS PRINCIPLES 


For the case of turbulent, viscous flow around the 
spherical bulbs shown in Figure (2), a complete solution of 
the Navier-Stokes equations is required to obtain accurate 
fluid velocities. Due to the difficulty in solving the 
complete Navier-Stokes equations, a commercially available 
computational fluid dynamics program has been employed to 
determine the velocities in a rectangular region around a 
Single bulb. Figure(3) shows a standard K-€ method solution 
of the Navier-Stokes equations for flow of air around a 
Circular cylinder generated by the FLUENT program. The 
shape of a 0.03m diameter circular cylinder in a 0.05m x 
0.20m flow field has been approximated by a series of grid 
points (100 wide by 25 high). The program calculates x and 
y velocities at every grid point in the flow for input into 
the convection terms of Equation [1]. 

The following parameters were input to the FLUENT 


program (properties of air/solution parameters): 


Properties of air only assumed (no air/form. mixture) 
Length of domain (x direction, m)- 0.2 

Height of domain (y direction, m)- 0.05 

Diameter of Bulb (m)- 0.03 


Gate Semew(<mGaiie sy dir)— 100 x 25 
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Figure (3) - FLUENT Generated Streamlines For Flow of 


Air Around a UV Lamp in a Uniform Flow Field 
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ile 
Reynolds Number, Re, - 4,500 
Inlet Temperature (°C)- 35 
Inlet Pressure (Kpa)- 101.33 
Inlet Density (kg/m) - ee 93 
imeem vmseocity (kq/m—s)— 1.72 x 107° 
Inlet Turbulence Intensity (%)- 10.0 


Inlet Velocity (m/s)- 2.0 


Since it is beyond the scope of this report to examine 
the detailed accuracy of the FLUENT program output, these 
results will be assumed to be sufficiently accurate for the 
purpose of making design decisions. There are a variety of 
models and techniques which can be employed to increase and 
verify the accuracy of the flow analysis. However, since it 
is the purpose of this report to model the photolytic 
reactor ina fashion which well-exceeds the accuracy of 
Simple plug flow models which can be found in the 


literature, the unmodified FLUENT output is used. 
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Chapter 5 


CHEMICAL PRINCIPLES 


The chemistry associated with the photolytic reactor 
under consideration here is much the same as that of 
atmospheric chemistry with the exception that energy is 
added with a mercury lamp which emits several discrete 
wavelengths (predominantly 253nm) versus sunlight which 
contains these and other wavelengths. The chemistry of the 
photolytic reactor can be broken down into two categories 


including photochemistry and free radical chemistry. 


5.1 PHOTOCHEMISTRY 
Photochemical processes initiated by UV radiation begin 
with the absorption of a photon by a molecule as follows: 


K 
A + hv - A* [oi 


where, A“ is an exited A molecule and hv is a photon of UV 
energy. Upon excitation the A* molecule can undergo the 


following reactions: 


K 
— i, 
Dissoectation: A*— Bie [6] 


(A* broken down into stable molecules by itself) 
Ko 
Direct Reaction: A*+B-moDteéE [7] 
(A* combines with B to form two stable molecules) 
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K 
5 
Fluorescence: A* — A + hw [8] 


(A* returns to orig. state after release of energy) 


Collisional Kg 
Deactivation: A* +BmoaA+t+B + hv [9] 


CeeGEeeibasS tO OL1G. 8State after collision with B) 


Each one of the above reactions has a probability of 
occurrence called a quantum yield, ®,. For example, a 
quantum yield of ®, = 0.5 implies that there is 50% chance 
that A* molecules will undergo dissociation (Equation [6]) 
while there is a 50% chance A* molecules will directly 
react, undergo fluorescence or deactivate from a collision. 
Since all the A* molecules must undergo one of the above 


reactions: 
=e [10] 


The production of A* molecules is given by the 


following rate equation: 


ad[A*] = K[{A] Po 
Gi 


In Equation [5], since one A* molecule consumes one 


molecule of A: 


adfA] = -K[A] {de 
at 





is 


The symbol K is the overall "specific absorption rate" 
constant (also given in Equation [2]) and [A*], [A] are the 
concentrations of A* and A, respectively. Note that the 
overall rate constant corresponds to all possible quantum 
yields for A*. By combining equations [5] through [9], the 
rate equations corresponding to the dissociation, direct 
reaction, fluorescence, and collisional deactivation 
reactions (termed the first order rate equations) can be 


written as follows: 


G{[A] = -D,K[A] = -K, [A] (reduction in A due to [ess | 
dt dissociation) 

d[A] = -O)K[A] = -K [A] (reduction in A due to [14] 
dt direct reaction) 

dqd{A] = —®3K [A] = -K3[A] (reduction in A due to [eS] 
dt fluorescence) 

d[A] = -O,K[(A] = -K,[A] (reduction in A due to [16] 
dt collisional deactivation) 


The quantities K, (=@,K), Ky K3, and Ky, inesecoations [13] 
through [16] are termed the first order rate constants. 
Note that multiplying both sides of Equation [10] by K 

a reomrietetek5 + Kat Kh, = K. Therefore, the overall 
reaction rate for the reduction of A is equal to the sum of 
the individual rates for reactions which can occur upon 
excitation of A molecules. 


For the general case of atmospheric chemistry, the 





iL, 
first order rate constants can be determined through the 


following equation: 


A 
K, = \ o(A,T)®, (A,T) I (A) dA [17] 
A, 
O(A,T) = absorption cross section of the A molecule at 


operating temp., pressure and UV wavelength 
(cm*) . 

®, (1,T) = quantum yield or probability that molecule A will 
dissociate upon excitation (recalling that one 
molecule of A produces one molecule of A*) (K5 


would correspond to direct reaction rate). 
I(A)dA = irradiance of incident UV light (photons/cm*-sec) . 


ies et hhou Order race COonStant, K, corresponds to 
the rate at which molecule A will undergo dissociation (one 
of the four reactions mentioned above) when exposed to a 
broad spectrum of light energy in the atmosphere. 

In order to calculate the first order rate constants 
for combinations of specific, discrete wavelengths the 


following approximation is utilized: 


n 
K, = Le o(hy, T)@, (Ag, T) I (Ag) AA [18] 
7] 





ZO 
The quantities 6 and @®, above are average values at A;. The 
quantity I(A,)AA is sometimes expressed as J(A,) (actinic 


flux) which results in the following: 


n 
Csi 


5.1.1 DETERMINATION OF FORMALDEHYDE PHOTOLYTIC RATE 


CONSTANTS 


Formaldehyde participates in the following photolytic 


reactions (A<300nm): 


K 
aL 

HCHO + hv — H + HCO (HCHO is formaldehyde) L210) 
Ro 


The wavelength of UV light utilized in this device is 


predominately A = 253 nm. The bulb dimensions are as 


follows: 
D (diameter) = 0.033 m 
L (lemgen) = 1.5536 m 
Ewepeeucoca. bulb power mequirement) = 40 W (J/s) 


Perea) (o> CUmeme 1m visible spectrum per 


manufacturer) = 40% Pi .+4) 


Poa. (bulb Oueout aeeZos nm) = 95% Peete 





Ze. 


As an example, the first order rate constant at A=253nm 
(corresponding to the reaction of equation [20] above) is 


determined as follows: 


PI 
I 


hc/A = 6. 6256x1 0" °*J-s/photon (2. 9979x10°m/s) 
253210 mM 


7.8509 x 107293/photon 


Ayes = BDL = (33mm) (1553. 6mm) = 1611 cmé 


Pos3 = 40W(40%) (95%)= 14.44 J/s 


J (A=253nm) Poe3 = 14.44 J/s 


Apyrp(E) 1611 cm*(7.8509 x 10719s/photon) 


1.1419 x 101° photons/cm’s 


O(253nmy25C) = 3.79 x 107*1cm*/molecule (R. Atkinson, 
et al. (1992) ) 


zoom, 25C) = 0.3 (R. Atkinson, et al. (1992)) 


Ky O(A, T)@, (A, T) F(A) 


3.79x10721cm2/molecule (0.3) (1.1419x10?®phot . /cm*s) 


1.30 x 107° (1/s) 


K = overall specific absorption rate (yields both 


reactions 1,2) 





Be 


K = Ky + Kp = O(A,T)®, (A,T)T(A) + (A, T) Bo (A, T) (A) 


O(A,T) (By (A, T) +o (A, T)) F(A) 


ee eee oe) — U-> (R, Atkinson, et al. (1992)) 


K = (3.79x1072tem?/molecule) (0.3+0.5) x 
(1.1419x10?®photons/cm?s) 


= 3.47 x 107° (1/s) (yields both reactions 1,2 
(cicmeZO yel 2. |) ) 


5.1.2 LIGHT INTENSITY DISTRIBUTION THROUGHOUT THE REACTOR 
The rate calculations above are based on the intensity 
at the bulb surface. This formulation is based on the Beer- 
Lambert law applied to weak absorbance conditions such as 
the atmosphere (Finlayson-Pitts, B. J., and Pitts, J. N., 
Jr. (1986)). Since the energy added via the bulb per unit 
area and corresponding rate constants vary indirectly with 
the radial distance from the bulb surface, Equation [2] must 


be modified as follows: 


ee Ry a) C Z| 


radial distance from bulb centerline. 


where, li 


wacias or bulb. 


a 
Il 


In addition, the light intensity in an actual 
photochemical reactor is also attenuated according to the 


well known Beer-Lambert law (Finlayson-Pitts, B. J., and 
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Piece. , Jr. {1O986))which is written as follows: 


I = Io exp{-ONP} (53) 


where, Io = monochromatic source intensity (photons/cm*-sec) . 
O = absorption cross section (cm*/molecule) . 
N = molecules/cm? (molar concentration). 


P = path length from the source (cm). 
For low concentrations and/or low absorbance, as is the 
case with this model, the effects of concentration on light 


intensity can be shown to be negligible. For this reactor: 


3.79 x 1074tcm*/molecule 


Go = 
N = 1.229 x 104° molecules/cm? (equivalent to 50 PPM) 
P = approx. 20 cm 


Substituting into Equation [23] yields: 


Ti S| ae x 10721¢m2/molecule) x 


HH 
i 


(1.229 x 104° molecules/cm?) (20 aul 


0.99991 Io 


In view of the above result it is clear that the weak 


absorbance assumption is a reasonable one. 
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In order to accurately model a multi-component air/VOC 
mixture where light absorbance must be considered, the 


following multi-reactant Beer-Lambert law is utilized: 


I = Io exp{-[0,N, ote OoN>d aM DS ce ]L} [24] 


etmgmenis celatatonship an intensity function, I = £(x,y) 
can be determined using the actual concentration of 
materials in the reactor which are capable of absorbing 
ilies Fritz 

The specific photochemistry of reactants which could be 
considered in this model are contained in Appendix B. 
However, for the purpose of this initial model only the 
photolysis of formaldehyde is under consideration. 
Therefore, the intensity will be considered to be a function 
of the radial distance from the bulb surface only due to the 
weak absorbance of formaldehyde which was demonstrated 


above. 


5.2 FREE RADICAL CHEMISTRY 

As can be seen in Appendix B there are a multitude of 
reactions which can be included in this model as there are 
numerous materials normally contained in an airstream which 
are capable of undergoing photolysis and reacting with 
formaldehyde. However, due to the great complexity of 


including all these reactions, this initial model will be 





Z5 
only concerned with the removal of formaldehyde through 
photolysis alone. Therefore the "generation" term of 
Equation [1] will be that given by Equation [20] with no 
other reactions included. 

Note that in making the above simplification, reactions 
with free radicals have been ignored. Therefore, the 
results of this investigation apply only to photolytic 
reactions. This analysis is a logical first step toward a 
more comprehensive analysis that will include free radical 
reactions. This analysis also provides valuable insight 
into the photolysis of formaldehyde in a non-uniform flow 


field. 
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Chapter 6 


CLOSED FORM SOLUTION FOR SIMPLIFIED FLOW FIELD 


In order to get a general idea of the size of the 
actual pollution control device required to reduce a given 
formaldehyde input concentration by a factor of two through 
photolysis, a simplified flow field model can be employed. 


Restating the mass conservation equation: 


-vyoVc + DV*c +Q = dCc/dt [25] 


The following simplifying assumptions will be made: 


A. Two-dimensional, inviscid, steady flow. 


(v = velocity in x direction = constant, 
velocity in y direction = Q) 


B. Constant radiance intensity throughout the reactor. 
C. Formaldehyde photolysis only. 


D. Diffusion negligible. 


With these assumptions, the above equation simplifies to: 


—-veVC + Q= 0 [26] 


Since the only reaction involves the photolysis of 


formaldehyde, substitution of Equation [2] into Equation 





vam | 


(26], performing the dot product and simplifying yields: 


v dC = -KC [27] 
ax 


Rearranging terms and solving for C(x) yields: 


Pe) Saye.) = —AX/V [28] 
ie eC.) =~ rear [29] 
C/o = GxXpi akeg v } [30] 


The values assumed are: 


C/G, = O58 
K = 3.47 x 107° (1/s) (as calculated) 


ele, S 


Substituting in values and solving for x yields: 


Gea 7 eee eo COrmes ponding to a 50% C reduction 


= 18.9 miles 


This result is a very pessimistic one which predicts 
poor performance for a reactor which relies on the 
photolysis of formaldehyde alone. However, in examining the 
general solution of the analytical model for different 


values of V and K, it appears that the analytical model can 





Formaldehyde Concentration (PPM) 


100 
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Formaldehyde Concentration vs. Distance 
Along an deal Two-dimensional Photolytic 
Reactor for Various Velocities and Rate Constants 


ko, V-5 tbs 
K=Ko, V=2 (U's 
K=10Ko, V=S fis 
K=1ORa, V=2 fs 
Ko = 5.47 x 10-5 (i/s) 





10600 20000 30000 40000 50000 
X Distance (ft.) 


Figure (4) - Formaldehyde Concentrations in an Idealized 


Two-dimensional Columnar Photo-reactor 
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be significantly improved (see Figure(4)) by employing a 
higher intensity mercury lamp (increasing K) and/or by 
decreasing the flow velocity. 

Increasing the intensity of the lamp is a viable 
option. There are currently bulbs available which are 
capable of over 1000 times the power of the mercury vapor 
lamps which are used in the photolytic reactor of Figure 
(a) . 

Since a decrease in velocity beyond some practical 
value will only lead to a decrease in the amount air 
treated, increasing the bulb intensity is clearly the best 


practical alternative. 
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Chapter 7 


DEVELOPMENT OF THE NUMERICAL MODEL 


The numerical photolytic reactor model under 
consideration predicts the concentrations in a rectangular 
area around a typical bulb in the reactor of Figure (2). 
Since very little reduction in formaldehyde is expected in a 
small 21.2 x 6.6 cm domain utilizing the actual rate 
constant calculated in chapter 5.1.1, the rate constant has 
been increased 10,000 times to an attainable 0.35 (1/s) to 
achieve some discernable results within two decimal places. 
Although this bulb intensity/reaction rate increase is 
presently impractical from a cost standpoint (very expensive 
bulbs are required) it greatly improves the analysis and 
interpretation of the numerical results for the small flow 
domain under consideration. The following data describe the 


parameters of the model: 


length of the domain in x direction = .2121 Mm 


~ 


number of nodes in the x direction = 100 


~ 


Salze On SIncremenk: Ss in the x direction = 0.002143 m 


length of the domain in the y direction = .0664 m 


™ 


number of nodes in the y direction = 31 


~ 


SS ee (2 Ss Ss - 


, size of increments in the y direction = 0.002143 m 
Uo, inlet velocity = 2.0 m/s (const. across entire 


area) 





eek 
Gayl owecehonumaleGenvyde concentration = 50.0 PPM ( " ) 
Bulb diameter = .03 m 


K, overall rate constant at bulb = .35 l/s 


D, diffusion coefficient = 0.00002 m?/s 


Pee owEmonmeonvergence criteria = 0.000001 PPM 


(maximum difference between successive concentration 
calculations at each point at a given time) 


At, time increment = 0.01 s 


The computer model uses a finite difference technique 
in conjunction with a Gauss-Seidel iterative method to 
calculate the concentration field throughout the reactor at 
selected time intervals. The mass conservation equation is 


restated as follows: 


-~vyeVc + DV*c +Q = daC/dt [31] 


Since the convection terms dominate Equation [31], an 
upwind difference method has been incorporated into the 
above equation to improve the accuracy of the results. 
Simply stated, the upwind technique allows the concentration 
at any point i,j to be calculated based on points upwind 
Since concentration changes are much more closely related to 


the fluid flow than they are to diffusion and photolytic 
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effects. Hence, differencing in the upwind direction is 
utilized for the convection terms. For the sake of 
simplicity, central differencing is used for the diffusion 
terms. The convection term of Equation [31] will now be 
developed in terms of finite differences. 

The equation which follows, Equation [32], is an upwind 
difference expression which is valid when both the fluid 
Veocteymectmenents, V, and v, are positive. 


y 


—vyoVC = Savery 10) ce ey ey 4) (Ci 474, 5-1) (S21 


a 
When both the fluid velocity components are negative: 
ee ere Ca 4-Cy 4) — (Vy /h) (Cy 54.-Ci, 4) [33] 


An expression derived from Equations [32] and [33] above 
which will always yield an upwind difference for any 
combination of velocity component directions can be written 
as follows: 

=veVC = -(1/2h)| Cee CeCe) + 
(lvy|+ vy) (Cy, 5-Cy, gaa) + (¥y7 [¥K}) (Coa, 37Ca, 3) + 


While the above equation may appear unnecessarily lengthy, 


= 





So 
it greatly increases the speed of numerical calculations as 
a decision does not have to be made (FORTRAN IF statement) 
at each point in order to apply the correct equations based 
on the velocity component directions. Without condition 
statements in the calculation of the convective terms, much 
of the program can be vectorized to greatly reduce the 
computation time. 

In terms of central differences, the diffusion term can 


be written: 


DVc = (D/n?)l Cia cca 


“Si oa Ci yaa) fo! 
The remaining terms are: 


Own = “KC ie) 


ag) 
= n ih 

The superscripts n and n-1l above refer to the present 
value and the value at the previous time, respectively. It 
is important to note that a backward difference scheme has 
been employed for the time dependant solution of the 
governing equation. Although all terms on the left side of 
the governing equation (Eqn. [31]) have been developed using 


concentration values at the present time, the superscript n 
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has been omitted for clarity. The backward difference 
method results in an implicit set of equations which can be 
solved by the Gauss-Seidel iterative method to find the 
concentration at each point in the reactor at a particular 
point (iteration) in time, n. This method may seem 
cumbersome in comparison to the forward difference method 
which can explicitly calculate the concentrations at a 
future point in time without iteration. However, the 
backward difference method has higher stability which allows 
a much larger time step to be used. 

SUsseaem@erng Heuaewons [33] through [37] into Equation 


fo) and “Ssolvang for Cig yields the following iterative 


equation: 
Ci al = 1 (1/2) Ao oe) 
Perley + hK + (h/At) 
C nese iin if OSes es Gtpmesea ct + (|v|- 72) Ge gas + 
elie ge + MD ynricy.y, 5h + cgay, seth + 
eae + C54 + (EYABC Gar" [38] 


The superscripts k and k+l in equation [38] represent 
concentration values calculated at the past and present 


iMe=amtons,, respectively. Note that the C; and Ci 


Saas = 
terms have the superscript k+l whereas the Ci+1,j and 
Ci 441 terms have the superscript k. This is due the 


iterative “Search pattern" which sweeps from j=1 to 
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N (y=0 to y=H) for each i until i=M (x=L) is reached. This 
particular sweep pattern in conjunction with the Gauss- 


Seidel method makes use of the latest values (k+1) for 


Ci-1,3 and Ci 4-1 in calculating the updated Gan value. 
The boundary conditions of the flow domain are set as 
follows: 


he WEileee= CC. — Const. = Co (1=1;183225) [39] 
B. Top and bottom surfaces - dC/dy = 0 (symmetry) or, 
Ci 541 = Ci, 5-1 (aay j=l, 25) [40] 
eee eet > — 9GG/ax, dc/dy = 0 (solid boundary) or, 
CG 


= iS (see program, Appendix A) [41] 


gis ioe, 
Ci 441 = Ci 45-1 (See program, Appendix A) [42] 
D. Exit - dC/dx = 0 (minimal change at exit) or, 


ce « = c 


The computer program included in Appendix A uses 
Equation [38] (including the relationships of Equations [39] 
through [43] when required) to calculate and update new 
concentration values at all grid points in the domain per 
each iteration. The iterations converge to the solution for 
a given time step when the present and previous 
concentration values for each point in the entire domain 

k+1 


differ by less than an epsilon value (i.e., Ci4 -c, 4 


See torrall points im the d@main) . 


Figure (5) is a graphical representation of the program 
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y(m) 
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0.02 


<I 
48.10 ue 3 
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oe UE ne in ~. === = 
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y(a) 
0 .0o4 


- 0.02 





Figure (5) - Formaldehyde Concentrations Within the 


Numerical Photo~reactor Model 





of 
output for selected values of time. The solution at t = 10 
sec is considered to be a steady state solution. This was 
determined by comparing the solutions at t = 10 sec and t = 
100 sec which did not differ to two decimal places for the 
Vis@eliagosiiy Or DOlNLS in the domain. In addition, to 
ensure that a time-accurate, steady state solution was 
obtained, the computer program was run at At=0.005 s (half 
the original time step) out to a 10 sec solution. The 
results of this run did not differ from the original 10 sec 
solution to over 4 decimal places for the entire domain. 

The numerical concentration predictions for the "start- 
up" of the reactor make sense from a practical standpoint. 
The smallest concentrations are seen at the bulb surface and 
wake regions due to the low fluid velocities there. In 
addition, these regions of lesser concentration migrate 
further downstream with increased time as they are swept 
into the bulb wake region. 

The steady state isopleths also appear reasonable in 
form. As expected, these concentration contours generally 
form the same shape as the streamlines depicted in the 
FLUENT output due to the dominance of the convection terms. 
However, the isopleth shapes differ from the streamlines due 
to diffusive transport and photolysis. In the next chapter 
it will be shown by comparison with the simplified flowfield 
model that the magnitude of formaldehyde reductions shown in 


Figure (5) are reasonable predictions. 
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Chapter 8 


ANALYSIS OF RESULTS 


In order to obtain a general performance check on the 
numerical solution, the simplified flow model developed in 
Chapter 5 can be employed by inputting the numerical model 
parameters discussed in Chapter 7 into Equation [30] as 


follows: 


Co = 50 (PPM) 


ho Ooo (17s) 


Ze Omen, S) 


V 


x 0.2121 (m) 


SUesctituering Values into equation [30] yields: 


@) 
i 


Co exp{-Kx/v} 


50.0PPM exp{-0.35s7+(0.2121m) /2.0m/s} 


48.17 PPM 


The above result compares favorably with the 49.37 PPM 
average steady state concentration leaving the numerical 
reactor model which was found through a mass average taken 


across the exit as follows: 


C = v.C, Lv, [45] 


exit, avg. at 
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In the next paragraph the analytical and numerical results 
are compared in terms of destruction efficiencies. 

An additional comparison between the numerical and 
simplified flow field models is obtained by ic roducing 
parameters which define the destructive efficiency of the 
photolytic reactor. Since it is the objective of the 
reactor model to eliminate formaldehyde, the efficiency of 


Eieomwentltemrmeactor section can be defined as follows: 


Bele (Coxit, avg. /Cinlet) [46] 
In addition to the above equation it is also useful to 
define an efficiency which is based on the bulb height 


alone. This term is defined as follows: 


Bap =1- (Coxit,avg./Cinlet) [47] 
Equation [47] calculates an efficiency based on a collimated 
region with a height equal to the radius of the bulb. In 
Equation [47] the term Cexit,avg./Cinlet is taken from the 
centerline of the reactor up to the bulb radius. Table-1l 
compares the efficiency of the steady state numerical 
Simulation of the photolytic reactor section with that of 


the simplified flow field model. 
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Simplified flow 
Field Model 
(PPM) 


Numerical Model 













(PPM) 







C (average across 
entire reactor 
section) 


C (average up to 
bulb radius) 


E (based on entire 
reactor cross 
section) 















(based on bulb 


Fip 3 
width) 





Table-1 Comparison of Numerical and Simplified Flow Models 


For this set of conditions it appears that the 
Simplified flow field model is more efficient than the 
numerical reactor model when the efficiency is based on the 
entire reactor cross section. However, in examining the 
efficiencies based on the bulb cross section, the numerical 
reactor model which incorporates the actual flow conditions 
is superior to the simplified flow field model in its 
ability to destroy the formaldehyde. This result leads one 
to believe that a reactor containing a combination of bulbs 
which allows the incoming flow to "see" bulbs across the 


entire reactor area will out-perform the idealized columnar 
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reactor. 

The advantage of using multiple bulbs is best 
illustrated by considering a single bulb reactor of infinite 
height. In such a reactor the UV bulb would have little 
effect on the VOC concentrations in the far stream. As 
shown in Figure (5), the bulb has the most influence 
(concentration reductions) in a columnar region which is one 
to two bulb diameters in height. In other words, the 
_ incoming flow which exists outside this columnar region 
tends to escape untreated. Therefore, an arrangement of 
bulbs which guarantees that all the incoming flow will 
encounter a bulb cross section (1e., a staggered stack of 
bulbs/columnar regions in parallel with no "escape" regions 
in between) would have a destructive efficiency of no less 
than that given by Equation [47]. 

In addition to the general advantages of a parallel 
bulb arrangement as described above, large-scale destructive 
efficiency improvements for both series and parallel 
arrangements are predicted by drawing an analogy between 
particle collection and VOC destruction in a photolytic 
reactor. For such an analogy the multi-bulb efficiency is 
given by the following (Heinsohn, 1991): 


E 1 - E,)N [48] 


overall = 1 7 ( 


where, N is the number of bulbs. 
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A twenty-five (N = 25) bulb arrangement with single bulb 
Seetctcictes a5 Given by equacdonm(4/) (E, = 5.26% in this 


case) would have an overall destructive efficiency as 


follows: 
ee ee = 1 = (1 - 0.0526)-° 
= 74.13 


Although the particle collector/photolytic reactor 
analogy drawn here predicts very optimistic results, it will 
be necessary to confirm this relationship with a multi-bulb 
numerical model and/or experimental testing to achieve a 


high degree of confidence. 
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Chapter 9 


CONCLUSIONS 


This analysis clearly demonstrates the influence of 
actual flow conditions on the effectiveness of a photolytic 
reactor. Despite the simplifying (but reasonable) 
assumptions made in development of the numerical model, it 
provides valuable insight into the inner workings of the 
actual photolytic reactor. The graphical results of Figure 
(5) show that the velocity gradients, recirculation and 
iene which oceur within the reactor contribute greatly to 
the overall reduction of reactants. Since recirculation and 
mixing greatly increase the reduction of formaldehyde in the 
reactor, one would expect that unsteadiness in the flow 
field would tend to improve the results (noting that 
unsteadiness exists in the actual flow field under 
consideration due to vortex shedding). In addition, one 
would expect that an optimum Reynolds number, R, could be 
found (input velocity) which would maximize the benefits of 
mixing while still providing an acceptable level of air 
treatment. Therefore, future numerical optimization studies 
Should include both an unsteady flow field analysis and a 
comparison of results for different Reynolds numbers. 

Although the simplified flow field model developed in 
this report can be used to obtain a rough estimate of 


reactor performance, it is unable to predict the effects of 
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the actual flow field which were found to be of great 
importance through an analysis of the numerical model 
results. This fact reinforces the need to use a numerical 
model in the optimization of a photolytic reactor design. 

Comparison of efficiencies between the simplified flow 
field and numerical models provides some highly instructive 
results. Surprisingly, though the light intensity in the 
numerical model is more than five times less than that of 
the simplified flow field model for a large portion of the 
domain (due to the radial light attenuation from the lamp 
surface), the numerical model provides a comparable 
formaldehyde reduction performance when the efficiency is 
based on the entire reactor cross section. In addition, it 
appears that increasing the number of bulbs will greatly 
improve overall reactor effectiveness as indicated by a 
comparison of efficiencies based on the bulb width. Due to 
the optimistic results obtained in comparing the numerical 
and simplified flow field models, the next logical step for 
future research in this area should include the modeling of 
larger flow domains which contain a larger quantity of 
bus . 

This undertaking also showed that additional factors 
can readily be taken into consideration in future reactor 
modeling. The computer program of Appendix A and techniques 
described in this analysis can be extended to include 


additional chemical reactions and other factors such as 
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Beer-Lambert light attenuation by reactants which will 
result in a numerical model of much greater accuracy. 
Although this work is of limited scope for practical 
reasons, expansion of this analysis including the use of a 
high-speed supercomputer could achieve an accurate numerical 
Simulation of the entire photolytic reactor of Figure (2). 

Despite the good correlation between the numerical and 

Simplified flow field models developed in this report, a 
high degree of confidence can only be achieved through 
laboratory experiments whose results can be compared to the 
numerical ones. Therefore, future work should include 
experimental verification. This will be especially 
important when additional reactants and other factors are 
introduced. In the future, the benefits of a closed form 
solution as experienced in this development (simplified flow 
field model) will not be possible when modeling a highly 


complex reactor. 
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50 


C SET UP VARIABLES 


MeL ECIT DOUBLE PRECISION (A-H, 0-2) 
DereNe LON Vx(100, 33) , VY (700, 33) , FMC (100, 33) , FMCOLD (100, 33) 
X0=4.286 

yO=0 .0 

RO=1.50 

M=100 

N=31 

H=0.002143 

EPSI=0.0000001 

DT=0.01 

ITER=1 

ITERT=1 

ITEST=1 

ITERMAX=200 

ITERTMAX=10 


C ENTER FORMALDEHYDE CONSTANTS 


FMRCO=0.35 
FMD=0 .00002 
FMC0=50 .0 


C READ IN FLUENT GENERATED VELOCITIES 


C SET 


150 


OPEN (UNIT=5, FILE=’ SINGLEBULBLIS.DAT’ , STATUS=’ OLD’ ) 
K=1 

KK=1 

KKK=10 

READ (5,5) ((VX(I,33-d) , I=KK, KKK) , J=1,N) 
K=K+1 

KKK=K*10 

KK= (K*10) -9 

IF (K.LE.10) GO TO 2 

K=1 

KK=1 

KKK=10 

READ (5,5) ((VY(I,33-d) , 1=KK, KKK) , J=1,N) 
K=K+1 

KKK=K*10 

KK= (K*10) -9 

IF(K.LE.10) GO TO 3 

CLOSE (UNIT=5) 

FORMAT (5X, 10E12.4) 

THE INITIAL FORMALDEHYDE CONCENTRATIONS 
DO 10 I=1,M 

DO 10 J=1,N+2 

FMC (I,J) =FMCO 

FMCOLD (I, J) =FMCO 

CONTINUE 


C CALCULATE THE CONCENTRATIONS AT BACH POINT USING GAUSS-SEIDEL METHOD 


I 


DO 50 I=2,M-1 
DOr 50° 0=2N+1 
TEMP=FMC (I,J) 


C CALCULATE THE PHOTOLYTIC RATE CONSTANT 


FMRC=FMRC0*RO/ ( (100. *H* (I-1) -X0) **2.+(100.*H* (J-2) ) **2.+ 


$0.00000000001) **0.5 





oe 


C APPLY BOUNDARY CONDITIONS (BC’S) AT THE BULB (J=2,1=14; J=4,I=15; 
@ J=6, I=16, J=7, I=17, J=8, I=19) 
IF (((J.EQ.2) .AND. (I.EQ.14)) .OR. ((J.EQ.4) .AND. (I.EQ.15)).OR. 
$((J.EQ.6) .AND. (I.EQ.16)).OR. ((J.EQ.7) .AND. (I.EQ.17)).OR. 
$((J.EQ.8) .AND. (I.EQ.19))) 
SFMC (I,J)=(1/ (4. *FMD/H**2.+FMRC+1./DT) ) * 
S ((FMD/H**2.)*(2.*FMC (I-1,Jd)+2.*FMC(1I,J+1))+FMCOLD (I,J) /DT) 
C APPLY BOUNDARY CONDITIONS (BC’S) AT THE BULB (J=2,1=28;J=4,1I=27; 
C Joel —2o 70-7, 0—250,0—6, 1—23) 
IF (((J.EQ.2) .AND. (I.EQ.28)) .OR. ((J.EQ.4) .AND. (I.EQ.27)) .OR. 
$((J.EQ.6) .AND. (I1.EQ.26)) .OR. ((J.EQ.7) .AND. (I.EQ.25)).OR. 
StIeso©.8) AND. (I.EO.23))) 
SFMC (I, J)=(1/ (4. *FMD/H**2.+FMRC+1./DT) )* 
S ((FMD/H**2.)*(2.* FMC (I+1,J)+2.*FMC(I,J+1) ) +FMCOLD (I,J) /DT) 
C SKIP AREAS INSIDE THE BULB 
Pras tem 2) aN. (J.LB.3)) .AND. ((I.GB.15) .AND. (I.LE.27))).OR. 
(ao 2Gs 4). 2ND ws .LE.5)) sAND.«((I.GE.26) .AND. (I.LE.26))).OR. 
$((J.EQ.6) .AND. ((I.GE.17) .AND. (I.LE.25))).OR. ((J.EQ.7) . AND 
St Gae le) weanm. (t.LE.24)) )wOR. ((J.B0.8) .AND. ((I.GE.20) .AND 
S(i.lm.22)))) GO TO 50 
C CALCULATION FOR GENERAL POINTS 
IF ((J.GE.9).OR. 
(((I.LE.13) .OR. (I.GE.29)) .AND. (J.EQ.2)).OR. 
meters) .Of(feGr,.28)) AND. ((J.GE.3) .AND. (J.LE.4))).OR. 
(((I.LE.15) .OR. (I1.GE.27)) .AND. ((J.GE.5) .AND. (J.LE.6))).OR. 
(((I.LE.16) .OR. (I1.GE.26)) .AND. (J.EQ.7)).OR. 
(((I.LE.18) .OR. (I1.GE.24) ) .AND. (J.EQ.8) )) 
FMC (I,J) =(1/ ( (ABS (Vx (I,J) ) +ABS (VY (I,J) )+4.*FMD/H) /H+FMRC+1./DT))* 
S(( (ABS (VX (I,J) )+VX(1I,J) ) *FMC (I-1, J) + (ABS (VX (I,J) ) -VX(I,J))* 
SFMC (I+1, J) + (ABS (VY (I, J) ) +V¥(I,d)) *PMC (I, J-1) 4 
S (ABS (VY (I,J) ) -VY (I, J) ) *FMC (I, J+1))/(2.*H) + (FMD/H**2.)* 
S$ (FMC (I-1, 7) +FMC (I+1, J) +FMC (I, J-1) +FMC (I, J+1) ) +FMCOLD (I, J) /DT) 
IF (ABS (FMC (I,J) -TEMP) .GT.EPSI) ITEST=ITEST*0. 
50 CONTINUE 
C ASSIGN VALUES REQUIRED DUE TO BOUNDARY CONDITIONS 
DO 52 I=2,M-1 
FMC (I,1)=FMC (I, 3) 
FMC (I, 33) =FMC (I, 31) 
52 CONTINUE 
FMC (20, 8) =FMC (20,10) 
FMC (21, 8) =FMC (21,10) 
FMC (22, 8) =FMC (22,10) 
FMC (18, 7) =FMC (18, 9) 
EMC (24,7) =FMC (24, 9) 
FMC (16,5) =FMC (14,5) 
FMC (26,5) =FMC (28,5) 
FMC (15,3) =FMC (13, 3) 
FMC (27,3) =FMC (29, 3) 
C APPLY BC’S AT END (DC/DX=0) 
DO 55 J=1,N+2 
55 FMC (M, J) =FMC (M-1,J) 
IF (ITER.GE.ITERMAX) GO TO 60 
IF (ITEST.EQ.1) GO TO 60 


$ 
$ 
$ 
$ 
$ 
$ 





60 


62 


52 
ITER=ITER+1 
ITEST=1 
GO TO 15 
IF (ITERT.EQ.ITERTMAX) GO TO 64 
DO 62 I=2,M 
DO 62 J=2,N+1 
FMCOLD (I, J) =FMC (I,J) 
ITER=1 
ITERT=ITERT+1 
ITEST=1 
GO TO 15 


C WRITE OUT RESULTS 


64 


OPEN (UNIT=6, FILE=’ NEWSINGLE. OUT’ , STATUS=’ NEW’ ) 
OPEN (UNIT=8,FILE=’ SINGLEPLOT.DAT’ , STATUS=’ NEW’ ) 
DOPL00 I=15,27 

DO 100 J=2,8 


C ASSIGN VALUES TO ACHIEVE A DECENT CONTOUR PLOT USING GNUPLOT 


100 


Pee ecm perp, (J.bB.3)) .AND. ((1.GH.15) .AND. (I.LE.27))).OR. 


pee geemee a sean. (J -LE 5) ) .AND. ((1.GE.16) .AND. (I.LE.26))) .OR. 
5 ((J.EQ.6) .AND. ((I.GE.17) .AND. (I.LE.25))).OR. ((J.EQ.7) . AND. 
5 ((I.GE.18) .AND. (I.LE.24))) .OR. ((J.EQ.8) .AND. ((I.GE.20) . AND. 
Sugmeren.22)))) HMC(ZI,J)=30. 


SUM=0. 


C CALCULATE EXIT AVERAGES 


80 


sal 


5D 


66 


67 
68 
Go 
70 
qe. 


SUMD=0. 

SUM1=0. 

SUM1D=0. 

DO 80 J=2,N+t1 

SUMD=SUMD+VX (100,J) 
SUM=SUM+FMC (100,J9) *VX (100,37) 
SUM=SUM/SUMD 

DO 81 J=2,9 

SUM1D=SUM1D+VxX (100,3J) 
SUM1=SUM1+FMC (100, J) *VX(100,J) 
SUM1=SUM1/SUM1D 

K=1 

KK=1 

KKK=10 

WRITE(6,71) ITER, ITERT, SUM, SUM1,DT,EPSI 

WRITE (6,70) ((FMC(I,33-J) , I=KK, KKK) ,J=1,N) 
=K+1 

KKK=K*10 

KK= (K*10) -9 

IF (K.LE.10) GO TO 65 

DO 67, J=2,Ntl1 

BO 66, I=1-M 

WRITE (8,68) FMC(I,J) 
WRITE (8, 69) 

CONTINUE 

FORMAT (F12. 8) 

FORMAT (1X) 

FORMAT (10F12.8) 
FORMAT (216, 4E12. 4) 

STOP 

END 
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WELL-MIXED MODEL OF 
PHOTO-OXIDATION AND RADICAL OXIDATION PROCESSES 


Contents: 

Photolysis Rates for 337 nm... De Wy eoncncecss0c ccscstsscseserseseccsees Z 
Pleteiysis Rates fer 253 69285 nim.......(UV bulb)....................cccsssssees 3 
Photolysis Rates for laser & UV bulb................... (GGMMDIO) ...5.0eccccse0s0:.- 4, 
SORT SI 0 22 ----2-<-0 onaceosncnnsa0no0c csusgisnoennensignageasenenon ancoeddluleceuegerencasencens 5 
Reactions: | 

Ba eI ad ee 22 a5, one cee 2c og SsUUAECE ES Ga. addsesunses-eesoseseceaseneesseoees 6 

Methanol 

Carbon Dioxide 

I tM ooo, <.asecuees-ooesscldlinyesc1=c0daesoesccececsesescosseecereccceaseccsoeees fi 

Ozone : 

Fee EON) Tyce ccs nnsc-ca-coscocecreosaesasiddsisstnre-as-esssorseaseadscaceesenssassassesss 8 

once oo eee s ccs Re EEE ose ncnccsceccsnsaicsesstecsececcocnceces 3, 

H202 

Eira re a oe Pan 2 osc ens. seccnesscvaceccacadesncsdcescsoreasees 10 

HONO 

Ny COP Sul IRZCU EZ UISSO) 5 | a ll 

ely ote) lene acy IAC AIS. 1) 2) .-...2<<ccauswsarunacs...0c0.sceeocccesooccesencccecsoeces 12 

Carbon Monoxide (CO) 

One ah ONC C00 Oe | 13 

Formyl Radical (HCO) | 
PUMA ISTIC GAGCS ooo ccc dcsctececc--2--cecncencoaccocerocscosessccencececccocecceccocerecsececsaccece 14 
oe cern 5 or 15 

Note: 


R represents a reaction; P represents a photolvtic reaction; KR 
represents a rate constant; KP represents a photolytic rate constant. 





Photolysis: 
x for the laser experiment will be 337 nm 
input file "RATE" 
Eppoton=590*10-2! J/Photon 


Actinic Flux=20*1015 Photons/(cm2*s) 


HCHO+hv — H*+HCO* 4<300nm 
HCHO+hv — CO+H? 42300nm 
CH30H+hv — negligible 

O3+hv — O(1D)+02 2s411nm 
O3+hv — O(3P)+02 4<1180nm 
O2+hv — O(1D)+0(3P) 2<175nm 
O2+hv — O(3P)+O0(3P) 2<242nm 
O2+hv ~ O(1D)+0(!D) A<137nm 
H202+hv — OH*+OH* 2<557nm 
H202+hv — H20+0(1D) <359nm 
H202+hv — H*+HO?2* 2<324nm 
H202+hv — 2H*+0? A<Z13nm 
NOz+hv — NO+0(1D) As244nm 
NOv+hv — NO+0(3P) 2<39§nm 
H2O0+hv — H2+0(3P) 4<243nm 
H20+hv —H*+OH* A<239nm 
H20+hv — H2+0(1D) A<176nm 
CO2+hv — CO+O0(1D) 

COQr+hv — CO+0(3P) 

HONO+hv — NO+OH* <591nm 
HONO+hv — H*+NO?2 4<367nm 
HONO+hv — HNO*+0(3P) 4<283nm 


4.08E-5 sl 
2.Oinieeaes- L 
O 


O 

Q, (? = QO) 
0 
0 


0 


1.40E-5 s-1 

0, (o = O) 

O, (> = O) 
0 


O 
7.08E-3 sl 


0 
0 
0 


negligible 
negligible 


3.23E-3s-1 


O, (& = O) 
0 


[P1A] 
(P1B] 
[P2] 


[P3A] 
[P3B} 
[P4A] 
(P4B] 


[P4C] 


[PSA] 
[P5B} 
[PSC] 
[PSD] 
[POA] 
[PGB} 
[PTA] 
[7B] 
[P7C] 


[P50A] 
[P50B} 
[PS50C] 
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Atkinson | 
Atkinson 


Atkinson 
Atkinson 
Atkinson 
Johnson/ 
Atkinson 

Johnson/ 
Atkinson 

Atkinson 
Atkinson 
Atkinson 
Atkinson 


Atkinson 
Atkinson — 


Atkinson 
Atkinson 
Atkinson 


Atkinson 
Atkinson 
Atkinson 





Photolysis: 


» for the UV Bulb experiment will be 253 nm and 185 nm 


input file "BULB" 
Ephoton=590%*10-4! J/Photon 


Actinic Flux: J(253)=1.14*1019 Photons/(cm2*s) 


J(185)=3.52*1014 Photons/(cm2+s) 


HCHO+hyv — H*+HCO* 
HCHO+hv — CO+H2 
CH30H+hv — negligible 


O3+hv — O(1D)+02 
O3+hv — O(3P)+02 
O2+hv — O(1D)+O(3P) 
O2+hv — O(3P)+O(3P) 


O2+hv + O(4D)+0(4D) 


H202+hv — OH*+OH* 
H207+hv — H20+0(!D) 
H202+hv — H*+HO?* 
H202+hv — 2H*+0O? 


NO2+hv — NO+O0(!D) 
NO2+hv — NO+0(3P) 


H20+hv —H2+0(P) 
H20+hv — H*+OH* 
H20+hv — H2+0(1D) 


COQ2+hv — CO+0(1D) 
CO7+hv —CO+0(3P) 


HONO-+hv — NO+OH* 
HONO+hv — H*+NQO?2 
HONO+hv — HNO*+0(3P) 


A<300nm 
A2300nm 


<41inm 
2<1180nm 
A<17Snm 

A<242nm 


A<137nm 


A<55/nm 
A<359nm 
A<324nm 
4<213nm 


Asl 44nm 
A<398nm 


A<?243nm 
A<?239nm 
A<176nm 


A<591nm 
A<367nm 
A<283nm 


1.30E-5s-! =[P1A] 
2.16E5s! [P1B] 
[P2] 
1.12E-1s1 [P3A] 
1.32E-2 sl — [P3B] 
0 [P4A] 
1.23E-9 st [P4B] 
0 [P4C] 
8.38E-4s71 [PSA] 
O, (> = 0) [P5B] 
O,(p=0) ~~ [P5C] 
O, (> = O) (P5D}] 
Ciie=0) "WRex| 
2.47E-4 el [PGB] 
O,(b=0) ([P7A] 
1.09E-5 st  [P7B] 
O [P7C] 
negligible 
negligible 
1.69E-3 s} [PSOA] 
O,(@=0) [PS5OB] 
0,(@=0) [P50C] 
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Atkinson 
Atkinson 


Atkinson 
Atkinson 
Atkinson 
Johnson/ 
Atkinson 
Johnson/ 
Athimson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 


Atkinson 
Atkinson 


Atkinson 
Atkinson 
Atkinson 


Atkinson 
Atkinson 
Atkinson 





Photolysis: 


a7 


» for the UV Bulb and the laser experiment will be 337, 253, and 185 nm 


input file "COMBO" 


Ephoton=>) 90*10-41 J7Photon 


Actinic Flux J(337)=20.0*1015 Photons/(cm2*s) 
J(253)=1.14*1016 Photons/(cm?*s) 
J(185)=3.52*1014 Photons/(cm2*s) 


HCHO+hv — H*+HCO* A<300nm 1.30E-5s-! [P1A] Atkinson 
HCHO+hv — CO+H? A2300nm 3.86E-5 sl} [PIB] Atkinson 
CH30H+hv — negligible fe =e - 
O3+hv ~ O(!1D)+02 As41linm 1.12F-1s°! [P3A] Atkinson 
O3+hv — O(3P)+02 A<1180nm 1.32E-2s°1 [P3B] Atkinson 
O2+hv + O(1D)+0(3P) .<175nm 0 [P4A] Atkinson 
O2+hv— O(3P)+0(3P) a<242nm —=1.23E-9 s*!_ [P4B] _ Johnson/ 
Atkinson 
O2+hv ~ O(4D)+0(1D) h<137nm O [P4C]  Johnson/ 
Atkinson 
H202+hv — OH*+OH* aA<557nm 8.52E-4s°1 [PSA] Atkinson 
H202+hv — H20+0(1D) 2<359nm O,(o=0) [PSB] Atkinson 
H202+hv — H*+HO?* 2<324nm O,(@=0) ([P5C] Atkinson 
H202+hv — 2H*+0O2 ~ 4<213nm O,(b=0) ([P5D] Atkinson 
NO2+hv — NO+O(1D) 4s244nm O,(o=0) [PGA] Atkinson 
NOz+hv — NO+0(3P) 4<398nm 7.42E-4s-1 [PGB] Atkinson 
H2O+hy —H>+0(3P) .<243nm 0,(6=0) ([P7A] Atkinson 
H20+hv — H*+OH* 2<239nm 1.09E-5 sl [P7B] Atkinson 
_ H2O+hv +H2+0(1D) A<176nm 0 [P7C] Atkinson 
CO2+hv > CO+0(!D) negligible 
CO2+hv — CO+0(3P) negligible 
HONO-+hv — NO+0OH* 4<59Inm 2.61E-3s°1 [P50A] Atkinson 
HONO+hv — H*+NO2 2<367nm O,(o=0) [PS50B] Atkinson 
HONO+hv — HNO*+0(3P) 2<283nm 0,(@=0) [P50C] Atkinson — 





OH*+03=HO2*+02 


HO2*+03=O0H*+202 


NO+03=NO2+02 
NOQ2+03=NO3+02 


HCO*+02=HO2*+C0 
HCHO+OH*=H*+CO+H20 
OH*+OH*=H20+0(1D) 
OH*+HO2*=H20+02 
OH*+H202=H20+HO2* 


NO2+OH*=HNO3 
CO+0OH*=C07+H* 
N+OH*=NO+H* 


HCHO+0OH*= H20+HCO 


H*+O2 = HO2* 
H*+O3 = OH* + O2 


HO2*+HO7*=H207+0? 
NO+HO7*=NO7+0OH* 
NO2+HO2*=HONO+0? 
CO+HO2*=CO2+0H* 
O(2D)+N2=O0(3P)4N2 
O(3P)+02+M=034+M 
NO+0(3P)+M=NO7+M 


NO27+0(3P)=07+NO 
O(3P)+03=02+02 


O(1D)+H20=OH*+0H* 
O(1D)+03=02+0(3P)+O(P) 
O(1D)+02=O0(5P)+02 


N+02=NO+0(3P) 
N+03=NO+0 
NO+OH*=HONO 


HONO+0H*=H20+NO? 
CH30H+0H*=0.85 {CH20H*}+ 
0.15{CH30*+H20} 
CH2OH*+02=HCHO+HO?2* 
CH30*+02=HCHO+HO>* 


6.7*10714 
2.0*10715 
ieseOs! 4 
1.2*10713 
5.6*10712 
1at0-!1 
1.9*10712 
1.1*10-10 
1.7*10°12 
1171107)! 
2.9*107!3 
4.9*10°11 
1.0*10°7!! 
1.2*10°12 
2.9*10°11 
1.6*10°!2 
8.3*10712 
1.4*10714 
1.9*10°32 
1.8*10°11 
2.9*10°11 
3.0*10-11 
9.7*10712 
9.5*10°15 
2.2*10719 
1.2*107}1 
4.0*107!! 
8.9*10717 
1.0*10716 
1.1*10°711 
4.9*1Q°12 
7.8*10713 


9.8*10712 
1.9*10715 


Radical Chemistry: 


cm3 molecule“! s-1 
cm> molecule! s-1 
cm3 molecule”! s-1 
cm? molecule™! s-1 
cm3 molecule”! s-1 
cm3 molecule’! s-1 
cm3 molecule’! s71 
cm3 molecule’! s-1 
cm3 molecule! s-1 
cm3 molecule”! s-1 
cm? molecule’! s-1 
cm> molecule! s-1 
cm3 molecule! s-1 
cm3 molecule’! s-1 
cm? molecule”! s-1 
cm3 molecule’! s-1 
cm? molecule"! s-1 
cm? molecule"! s-1 
cm? molecule’! s-1 
cm> molecule’! s-1 
cm3 molecule’! s-1 
cm3 molecule’! s-1 
cm3 molecule! s-1 
cm3 molecule’! s-1 
cm? molecule”! s-1 
cm3 molecule-! s-1 
cm? molecule’! s-1 


cm3 molecule’! s-1 


cm? molecule-! s-1 
cm? molecule-! s-1 
cm3 molecule7! s-1 
cm? molecule’! s-! 


cm? molecule! s-1 
cm3 molecule”! s-1 
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DeMore 
DeMore 
DeMore 
DeMore 
Anatasi 
Seinfeld 
DeMore 
Atkinson 
DeMore 
Seinfeld 
DeMore 
Atkinson 
Atkinson 
Anatasi 
Anatasi 
DeMore 
DeMore 
Stockwell 
Atkinson 
DeMore 
Seinfeld 
Atkinson 
DeMore 
DeMore 
DeMore 
DeMore 
DeMore 
Atkinson 
Atkinson 
Atkinson 
Atkinson 


Atkinson 


Atkinson 
Atkinson 





Formaldehyde Reactions (HCHO): 


Photolysis: 


HCHO+hyv —--H*+HCO* -4<300nm 1.30E-5s 7} 
HCHO+hv --CO+H2 A=300nm 3.86E5 si 
Radical Chemistrv: 

HCHO+OH*=H*+CO+H20 1.1*10-1! cm3 molecule! s-1 
HCHO+OH*=H20+HCO* 1.0*10-!! cm3 molecule7! s-1 
CH20H*+07=HCHO+HO2* 9.8*10°12 cm3 molecule! s-1 
CH30*+02=HCHO+HO?2* 1.9*10°15 cm3 molecule! s-1 


d{[HCHO]/dt=-P1A-P1B-R10-R17+R61+R62 


Methanol Reactions (CH30H): 


Photolvsis: 
CH30H+hv — negligible 


Radical Chemistrv: 


CH30H+0OH*=0.85 {CH20H*}+ 7.8*10713 cm3 molecule’! s-1 
0.15{CH30*+H20} 


d{(CH30H]/dt=R6O 


Carbon Dioxide Reactions (CO2): 
Radical Chemistry: 


CO + OH*— CO? + H* 2.2*10-13 cm3 molecule7! s-1 
CO + HO2*-> CO? + OH* 1.9*10°32 cm3 molecule’! s-! 


d[CO2]/dt=P15+P24 


[P1A] 
[P1B] 


[R10} 
[R17] 
[R61] 
[R62] 


[P2] 


[ROO] 


[R15] 
[R24] 


oP, 


Seinfeld 
Seinfeld 


Seinfeld 
Atkinson 
Atkinson 
Atkinson 


Atkinson 


DeMore 
Atkinson 





Atomic Oxygen Reactions (QO): 


Photolvsis: 

O3+hv -O(1D)+02 4<411nm 
O3+hv —-O(3P)+02 a<1180nm 
O2+hv -O(1D)+0(3P) 4<175nm 

O2+hv -O(3P)+0(3P) 1<242nm 

On+hv -O(!D)+0(1D) 4<137nm 

NO2+hv -NO+0(!D) 4<244nm 
NO2?+hv —NO+0(3P) 1<398nm 

H2O0+hv —» H7+0(3P) A<243nm 

H20+hv — H2+0(!D) A<176nm 

H207+hv — H20+0(1D) 4<359nm 

HONO+hv — HNO*+0(3P) 4<283nm 

Radical Chemistry: 

OH*+OH*=H20+0(!D) 1.9*10712 


O(1D)+N2=O0(3P)+N2 1.8*107!! 
O(3P)+02+M=034+M - 2.9*1Q711 
NO+0(3P)+M=NOQ7-M 3.0*10711 
NO27+0(3P)=02+NO 9.7*10712 
O(3P)+03=207 9.5*10715 
O(!D)+H20=20H* 2.2*10-10 
OQ(1D)+03=07+20(3P) peste! 
O(1D)+07=0(3P)+0? 4.0*10711 
N+027=NO+0(3P) 8.9*10°17 


a4Assumed the products of R11 were negligible. 


Use Steady State Assumption: - 


live) 57} 
1.32E-2 s-l 
0 
1.23E-9 sl 
0 
O, (> = O) 


7.42E-4 s71 


O, (> = 0) 
O 

O, ( = 0) 

O,(@ =0) 


cm3 molecule”! s-1 4 
cm? molecule~! s-1 
cm 3 molecule-! s-1 
cm? molecule"! s-1 
cm? molecule7! s-1 
cm? molecule™! s-1 
cm3 molecule! s-1 
cm? molecule! s-1 
cm? molecule7! s-1 
cm molecule! s-1 


[P3A] 
[P3B] 
[P4A] 
[P4B] 

[P4C] 

[PGA] 
[PGB] 
[P7A] 
[P7C] 

[P5B] 
[P50C] 


[R11] 
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Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 
Atkinson 


DeMore 
DeMore 
Seinfeld 
Atkinson 
DeMore 
DeMore 
DeMore 
DeMore 
DeMore 
Atkinson 


([0(3P)]ss=(P3B+P4A+P4B+PGB+P7A+P50C+R30+R37+R38+R39)/(KR3 st 


NO2*KR35*03+KR32*NO+KR31*O0 2) 


[0(1D)]ss=(P3A+P4A+P4C+PGA+P7 C+PSB+R] 1)/(KR30*N2+KR36*H20+ 


KR37*03+KR38*O72) 





NO Reactions: 


Photolvsis: 

NOz+hv -»>NO+0(1D) 
NQ2+hv —-NO+0(3P) 
HONO+hv —NO+0OH* 


Radical Chemistrv: 

NO + 03 = NO? + O2 

N + OH* = NO + H* 

NO + HO2* = NO? + OH* 
NO + O(3P) + M=NO?7 +M 
NO? + O(3P) = O02 + NO 
N+07=NO+0(3P) 
N+03=NO0+0? 
NO+OH*=HONO 


1.8*10-!4 
4.9*10°11 
8.3*10-12 
30a)! 
9.7*10°12 
8.9*10717 
1.0*10716 
1.1*10°11 


As244nm 
14<398nm 
4<591nm 


O, (> = 0) 
7.42E-4 sl 
2.61E-3 s71 


cm3 molecule’! s-1 
cm3 molecule"! s-1 
cm3 molecule’! s-1 
cm3 molecule’! s-1 
cm3 molecule’! s-1 
cm3 molecule”! s-1 
cm3 molecule”! s-1 
cm3 molecule’! s-1 


[PGA] 
[POB] 


[PSOA] 


[R3] 
[R16] 
[R22] 
[R32] 
[R33] 
[R39] 
[R40] 
[R51] 


d{NO]/dt=P6A+P6B+P50A-R3+R16-R22-R32+R33+R39+R40-R51 


Hydrogen Peroxide Reactions (H202): 


Photolvsis: 

H202? + hv— OH* + OH* 
H202 + hv—-H20 + O(1D) 
H202 + hv— H*+HO?* 
H202 + hv— 2H*+0? 


Radical Chemistrv: 
OH* + H202 = H20 + HO?2* 
HO 2* + HO2* = H202 +02 


A<5/77 nm 
~<359 nm 
A<324nm 
Ree yS pin 


d{H202]/dt=-PS5A-P5B-P5C-P5D-R13+R20 


1.7*10°12 cm3 molecule! s71 
1.6*10°12 cm3 molecule’! s-1 


[P5A] 
[P5B] 
[PSC] 
[PSD} 


[R13] 
[R20] 
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Atkinson 
Atkinson 
Atkinson 


DeMore 
Atkinson 
DeMore 
Atkinson 
DeMore 
Atkinson 
Atkinson 
Atkinson 


Atkinson 
Atkinson 
Atkinson 
Atkinson 


DeMore 
DeMore 





NO?2 Reactions: 


Photolysis: 

NOvz+hv -NO+0(!D) 
NQ?+hv —--NO+0(3P) 
HONO+hv—H*+NQ? 


Radical Chemistry: 
NO+03=NO2+02 
NOQ2+03=NO3+02 
NO2+OH*=HNO3 
NO+HO2*=NO7+OH* 
NO2+HO2*=HONO+02 
NO+0(3P)+M=NO7+M 
NO7+0(3P)=02+NO 
HONO+0H*=H20+NO? 


A1<s244nm O, (> = O) 
1<398nm 7.42E-4 s-1 
A4<367nm O, (> = 0) 


1.8*10°14 cm3 molecule’! s-1 
1.2*10713 cm3 molecule’! s-1 
1.1*107!! cm3 molecule’! s-! 
8.3*10°12 cm3 molecule’! s-1 
1.4*10°12 cm3 molecule! s-1 
3.0*10-!1 cm3 molecule”! s-! 
9.7*10°12 cm3 molecule”! s-! 
4.9*10712 cm3 molecule’! s-1 


[PGA] 
[POB] 
[PS OB] 


[R3] 
[R4] 
[R14] 
[R22] 
[R23] 
[R32] 
[R33] 
[R52] 


d{NO2]/dt=-PGA-P6B+P50B+R3-R4+-R14+R22-R23+R32-R33+R52 


Nitrous Acid Reactions (HONO): 


Photolvsis: 
HONO+hv —-NO+OH* ~ 
HONO+hyv—H*+NO? 


HONO+hy—-HNO*+2(5P) 


Radical Chemistrv: 
NO2+HO2*=HONO+02 
NO+OH*=HONO 
HONO-+-OH*=H20+NO2 


A<591nm 2.61E-3 sl 
<367nm O, (> = O) 
2<283nm O, (> = O) 


1.4*10°!4 cm3 molecule’! s71 
1.1*107!1 em3 molecule7! s-1 
4.9*10-12 em3 molecule’! s-1 


d{HONO]/dt=-P50A-P50B-P50C+R23+R51-R52 


[P50A] 
[P5 OB] 
[P50C] 


[R23] 
[R51] 
[R52] 
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Atkinson 
Atkinson 
Atkinson 


DeMore 
DeMore 
Seinfeld 
DeMore 
Stockwell 
Atkinson 
DeMore 
Atkinson 


Atkinson 
Atkinson 
Atkinson 


Stockwell 
Atkinson 
Atkinson 





Hydroxyl Radical Reactions (OH*): 


Photolvsis: 
H202+hv —OH*+OH* 
H20+hyv—-H*+OH* 
HONO+hyv —-NO+OH* 


Radical Chemistry: 
OH*-+03=HO2*+02 
HO2*+03=OH*+202 
HCHO+OH*=H*+CO+H20 
OH*+OH*=H20+0(1D) 
OH*+HO27*=H20+02 
OH*+H202=H20+HO?2* 
NO2+0OH*=HNO3 
CO+OH*=CO7+H* 
N+OH*=NO+H* 

H*+O3 = OH*+02 
NO+HO27*=NO7+0H* 
CO+HO2*=CO2+OH* 
O(1D)+H20=OH*+0H* 
NO+OH*=HONO 
HONO+OH*=H20+NO? 


CH30H+OH*=0.85{CH20H*}+ 
0.15{CH30*+H20} 


DAssumed [OH*] to be low, considered R11 to be negligible. 


4<55/nm 
A<239nm 
1<591nm 


8.52E-4 $71 
1.09E-5 s71 
2.61E-3 s7l 


6.7*10714 
2.0*10°715 
1 he Coe 
1.9*10712. 
1.1*10-19 
1.7*10712 
1.1*1071! 
22163) 
4.9*10-11 
2.9*10711 
8.3*10712 
1.9*10732 
2.2*10710 
1.1*10-!2 
6.6*10°12 
7.8*10713 


Use Steady State Assumption: 


[OH*]}ss=(P5A+P7B+PS50A+R2+R19+R22+R24+2*R36)/(KR1*O3+KR10* 


cm> molecule"! s-1 
cm3 molecule”! s-1 
cm? molecule’! s-1 
cm? molecule"! s-1 b 
cm? molecule”! s71 
cm} molecule”! s-1 
cm3 molecule"! s-! 
cm3 molecule7! s-1 
cm? molecule7! s-1 
cm? molecule”! s-1 
cm? molecule! s-1 
cm3 molecule’! s-1 
cm3 molecule”! s-1 
cm? molecule”! s-1 
cm3 molecule7! s-1 
cm3 molecule”! s-1 


[PSA] 
[P7B] 


[PSOA] 


[R1] 
[R2] 
[R10] 
[R11] 
[R12] 
[R13] 
[R14] 
[R15] 
[R16] 
[R19] 
[R22] 
[R24] 
[R36] 
[R51] 
[R52] 
[RGO} 


Atkinson 
Atkinson 
Atkinson 


DeMore 
DeMore 
Seinfeld 
DeMore 
Atkinson 
DeMore 
Seinfeld 
DeMore 
Atkinson 
Anatasi 
DeMore 
Atkinson 
DeMore 
Atkinson 
Atkinson 
Atkinson 


HCHO+KR12*HO2+KR13*H202+KR14*NO2+KR15*CO+KR16*N+ 
KR5 1*NO+KRS5 2*HONO+KRGO*CH 30H) 





Hydroperoxy! Radical Reactions (HO2*): 


Photolvsis: 
H202+hv—-H*+HO2* 


Radical Chemistrv: 
OH*+03=HO2*+0?2 
HO2*+03=OH*+20?2 
HCO*+07=HO?*+CO 
OH*+H202=H20+HO?* 
H*+0O2=HO?* 
HO2*+HO2*=H202+02 
NO+HO7*=NO7+OH* 
NO2+HO7*=HONO+0? 
CO+HO27*=CO2+0H* 


CAssumed [HO?2] to be low, considered R20 to be negligible. 


A<324nm 


6.7*10-14 
2.0*10°15 
5.6*10°12 
1.7*10°12 
1.2*10712 
1.6*10°12 
8.3*10°12 
1.4*10°14 
1.9*10°32 


Use Steady State Assumption: 


QO, (=0) 


cm3 molecule’! s-1 
cm3 molecule"! s-1 
cm3 molecule’! s-1 
cm3 molecule7! s-1 
cm3 molecule71 s-1 
cm> molecule”! s-1 ¢ 
cm3 molecule’! s-1 
cm3 molecule7! s-1 
cm3 molecule’! s-1 


[P5C] 


[R1] 
[R2] 
[R9] 
[R13] 
[R18] 
[R20] 
[R22] 
[R23] 
[R24] 
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Atkinson 


DeMore 
DeMore 
Anatasi 
Seinfeld 
Anatasi 
Seinfeld 
Seinfeld 
Stockwell 
Atkinson 


[HO2*]}ss=(P5C+R1+R9+R13+R18)/(KR2*03+KR22*NO+KR23*NO2+KR24*CO) 


Carbon Monoxide Reactions (CO): 


Photolvsis: 
HCHO+hv — CO+H? 


Radical Chemistrv: 
HCO*+07=HO?*+CO 


HCHO+OH*= H*+CO+H20 


CO+OH*= CO2+H* 
CO+HO?*= CO27+OH* 


4>300nm 


d[CO]/dt=P1B+R9+R10-R15-R24 


3.8G6E-5 s7l 


5.6*10-12 cm3 molecule7! s-1 
1.1*10°11 cm3 molecule’! s-1 
2.2*10°13 cm3 molecule’! s-1 
1.9*10°32 cm3 molecule"! sl 


[P1B] 


[R9] 
[R10] 
[R15] 
[R24] 


Atkinson 


Anatasl 
Seinfeld 
DeMore 


Atkinson 





Atomic Hydrogen Reactions (H*): 


Photolvsis: 
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HCHO+hv — H*+HCO* 4<300nm 1.30E5s-1 [P1A] Atkinson 
H207+hv —H*+HO2* 4<324nm O,(@=0) [P5C] Atkinson 
H202+hv — 2H*+0O2 4<213nm O,(@’=0) [PSD] Atkinson 
H20+hv — H*+OH* 4<239nm 1.09E-5 sl [P7B]) Atkinson 
HONO+hv —> H*+NO?2 4<367nm O,(@=0) [PS50OB] Atkinson 
Radical Chemistry: 

HCHO+OH*=H*+CO+H20 1.1*10-1! cm3 molecule7! s-1 [R10] Seinfeld 
CO+OH*=CO?+H* 2.9*10°!13 cm3 molecule! s-1 [R15] | DeMore 

N+OH*=NO+H* 4.9*10-!11 cm3 molecule"! s-1 [R16] Atkinson 
H*+O2 = HO?* 1.2*10°12 cm3 molecule7! s-1 [R18] Anatasi 

H*+O3 = OH*+ O2 2.9*10-!11 cm3 molecule”! s-1 [R19] Anatasi 

d{H*]/dt=P1A+P5C+P5B+P7B+P50B+R10+R15+R16-R18-R19 

Formyl Radical Reactions (HCO*): 

Photolvsis: 

HCHO+hyv — H*+HCO* 4<300nm 1.30E-5s-1 [PIA] Atkinson 
Radical Chemistrv: 

HCO*+02 = HO2*+CO 5.6*10-12 cm3 molecu 27! $71 [RI] Anatasi 


1.0*10-!1 cm3 molecule’! s-1 [R17] Atkinson 


HCHO+OH*= H20+HCO* 


d(HCO*]/dt=P1A-R9+R17 
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Atmospheric Gases: (Seinfeld, p. 8] 


Gas Avg. Concentration[ppm 

Ar 9340 

Ne 18 

Kr ial 

Xe 0.09 

N2 780,840 

OQ? 209,460 
CH4 1.65 

CO? 332 

© 0.05-0.2 

H2 0.58 

N20 0.33 

SO? 10-5-10-4 
NH3 10-4-10-3 

NO+NOQ?2 10-6-10-2 

O3 10-2-10-1 
HNO3 10-5-10-3 
OH 0.406*10-> ppb 
HO? 0.406*10-2 ppb 

O3 30 ppb 
A202 10 ppb 
H20 22,550 


He Syy 
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